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Abstract 

This paper formulates and studies a stochastic maximum principle for forward- 
backward stochastic Volterra integral equations (FBSVIEs in short), while the 
control area is assumed to be convex. Then a linear quadratic (LQ in short) 
problem for backward stochastic Volterra integral equations (BSVIEs in short) is 
present to illustrate the aforementioned optimal control problem. Motivated by the 
technical skills in solving above problem, a more convenient and briefer method 
for the unique solvability of M-solution for BSVIEs is proposed. At last, we will 
investigate a risk minimization problem by means of the maximum principle for 
FBSVIEs. Closed-form optimal portfolio is obtained in some special cases. 

Keywords: Forward-backward stochastic Volterra integral equations. Adapted 
M-solution, Optimal control. Stochastic maximum principle. Backward linear quadratic. 
Risk minimization problem 



1 Introduction 

Throughout this paper we assume that all uncertainties come from a common com- 
plete probability space (O, P) on which is defined a d-dimensional Wiener process 
(V7t)fg[o,T]- The main objective of this paper is to study the optimal control problem 
for the following forward-backward stochastic Volterra integral equation (FBSVIE in 
short) 

X{t) = ip{t)+ [ b{t,s,X{s))ds+ I a{t,s,X{s))dW{s), 

T ° ° T (1) 

Y{t) = ij{t)+ [ git,s,X{s),Yis),Zis,t))ds- [ Z{t, s)dW{s), 
Jt Jt 
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which for instance generahze the optimal control problems for stochastic Volterra inte- 
gral equations in [26] . 

The notion of M-solution for backward stochastic Volterra integral equation (BSVIE 
in short) of the form 



with t G [0,?^], was introduced by Yong in [20], which plays an important role in 
optimal control problem for stochastic Volterra integral equations (SVIEs for short). 
We refer the author to Lin [7], Yong ([M], [25]), Wang and Zhang ([21]) for a study 
of the wellposedness of BSVIEs in finite space, while Anh and Yong [1], Ren [18] in 
infinite space counterpart. 

One main feature of equation (2) lies in the dependence of the generator g on Z(s, t), 
and hence it is quite different from, more precisely, a natural generalization of the one in 
[7] and [2T]. Of course, the appearance of such term in g is not just means an extension 
from mathematical point of view, but also be of great importance in applications, (see 
Proposition 3.5 in [25] and Theorem 5.1 in [26].) It is interesting to realize that, so far 
as we know, it is just the term Z{s, t) rather than Z{t, s) in the generator g that plays 
a key role in both optimal control problem in [26] and dynamic risk measure in |25j . 

Optimal control of forward stochastic differential systems is a classical problem. 
When we consider the Pontryagin maximum principle for optimal controls of stochastic 
differential equations, the adjoint equation for variational state equation actually is a 
linear backward stochastic differential equation (BSDE for short). The wellposedness 
for nonlinear BSDEs was firstly studied by Pardoux and Peng |14j . Readers interested 
in an in-depth analysis of BSDEs can see the books of Ma and Yong [8], Yong and 
Zhou [28j and the survey paper of EI Karoui, Peng and Quenez [3]. As to the optimal 
control for stochastic differential equation, we refer the reader to, for example, Peng 
|15| for the general case of control domain being non-convex, and Yong and Zhou 
|28j for systematical analysis. On the other hand, optimal control for deterministic 
Volterra integral equation, particularly linear quadratic problem, was firstly studied by 
Vinokurov |19j . From then on some other extensions were developed, see, for example, 
[2], [3], [IT], [27] and the references cited therein. As to the stochastic version, Yong 
( |24j and |26j ) presented a maximum principle for SVIEs by means of BSVIEs, while the 
control is assumed to be convex. We also would like to mention the work of Oksendal 
and Zhang [13] in partial information setting without the help of BSVIEs. Along this, 
we will investigate the FBSVIEs case in this paper. To the best of our knowledge, so far 
little is known about maximum principle for FBSVIEs, and one aim of this manuscript 
is to close the gap. 

The scheme is designed around the three steps for FBSDEs in Peng [16], namely 
listing out the variational equation, obtaining the variational inequality and utilizing 
some key mathematical tools to finish the procedure. As we know, within the context 
of stochastic differential systems, Ito formula has received most attention largely due to 
its ad hoc role in many complicated calculations and proofs. For example one usually 




Y{t) = m + 



git,s,Yis),Zis,t),Z{t,s))ds- Z{t,s)dWis) (2) 



Jt 
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makes use of Ito formula in obtaining the convergence property for Xp and Yp (defined 
blew) for differential systems, see Lemma 4.1 in |16j . In fact, one key tool in deriving 
maximum principle in Peng [16] is just Ito formula too. Unfortunately, this efficient 
tool is failure in the Volterra integral systems and some related well properties are 
absent in this case. 

In this paper, new approaches are proposed to handle with the difficulties encoun- 
tered in above procedure. On the one hand, we will make use of the dual princi- 
ple, established by Yong in [26], for linear stochastic integral equation and its adjoint 
equation. Consequently we have to tackle four equations which perhaps means more 
mathematical expressions and notations involved after introducing another two more 
adjoint equations for FBSVIEs. As a result, it is our hope to choose appropriate form 
of adjoint equations so as to make the procedure as brief as possible. Fortunately, such 
adjoint equations really exist, see (10) and (22). On the other hand, we introduce a new 
equivalent norm for elements in '?^^[0,T], see (6), and use some common calculations 
and tricks employed in the conventional BSDEs case, thereby obtain some convergence 
results, which play a chief role in deducing the variational inequality. Notice that Ito 
formula does not appear in the above two aspects. 

Motivated by the new norm aforementioned, in the following we will provide a 
new method for the unique solvability of M-solution, which seems more convenient 
than the one in [26]. By the four steps in Theorem 3.7 in [26] we can see the process of 
constructing the M-solution clearly. From mathematical view, however, the whole proof 
is too complicated and uneasy to understand, which prompts us to seek an alternative 
one. We will carry out this course in detail in Section 3. 

A class of continuous time dynamic convex and coherent risk measures, perhaps 
allowing time-inconsistent preference unlike the conventional case, were introduced by 
Yong in [25] via BSVIEs of the form 

Y{t) = -m + r git, s, Y{s), Z{s, t))ds - C Z{t, s)dW{s), (3) 
Jt Jt 

with t G [0, T]. In the classical case, the terminal condition is usually a bounded random 
variable, representing the financial position at time T. However, in the situation under 
our consideration, we prefer to choose a process "0) representing the total wealth of 
certain portfolio process at time t which might be a combination of certain contingent 
claims, positions of stocks, mutual funds and bonds. Usually the '0 could be i3[0,T] (g) 
J^T-measurable, see the example on p. 13 in [25]. If we define a map q from L% [0, T] to 
L|[0,T] by Q{t]^l}{-)) = Y{t), with Y being the M-solution of BSVIE (3), given certain 
assumptions on it is shown in |25| that q could be a dynamic convex/coherent risk 
measure. The question is how to look for a appropriate portfolio that minimizes the 
risk of the wealth process ^ by means of the representation above in finance, i.e., to 
seek an optimal solution for the so-called risk minimization problem, see Mataramvura 
and Oksendal [9], Oksendal and Sulem |12j for more information on the above financial 
problem. We conclude this paper by giving a study of this problem with the help of 
the maximum principle. In some cases, the closed form of optimal portfolio is derived. 
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The remainder of paper is organized as follows. In Section 2, we give some prelim- 
inary results and notations which are needed in the following sections. A new method 
for the solvability of M-solution is presented in Section 3. We give the stochastic max- 
imum principle for FBSVIEs (1) as well as a backward linear quadratic problem in 
Section 4. At last, we investigate a risk minimization problem by means of maximum 
principle in the previous. Some explicit form solutions are derived. 

2 Preliminaries 

In this section, we will make some preliminaries. Let us specify some notation in 
this paper. For any R,S e [0,T], we denote A^[0,T] = € [0,T]2;i < s}, 

A[0,r] = G [0,r]2;t > s}. In what follows some spaces will be frequently used. 

Let L^^[0,T] be the set of the B{[0,T]) (g) Tt processes X : {0,T] x Q ^ R"^ such 

that E J^\X{t)\Pdt < oo. L^^[0,T] is the set of ah adapted processes X : [0,T] x 
^ ^ such that E \X{s)\Pds < oo. LP{0,T; L^O,T]) is the set of ah processes 
Z : [0,r]2 X $7 ^ i?"^'^' such that for almost all t G [0,r], Z{t,-) is /"-progressively 

P 

measurable satisfying E ( Jq \Z{t, s)\ ds j dt < oo. For notational clarity, we denote 

•HP[0, T] = ^^,[0, T] X LP(0, T- L%[i), T]). Next we shah cite the definition of M-solution 
introduced in [26] . 

Definition 2.1 A pair of (Y {■) , Z {■ , ■)) G ?^^[0,T] is called an adapted M-solution of 
BSVIE (2) on [0,T] if (2) holds in the usual Ro's sense for almost all t G [0,T] and, 
in addition, Y{t) = EY{t) + jl Z{t, s)dW {s) with t G [0,r]. 

The next two definitions are introduced by Yong in |25] . 

Definition 2.2 A mapping p : Ljr^[0, T] L|[0,T] is called a dynamic risk measure 
if the following hold: 

1) (Past independence) For any ^{■), "^(•) G L]r^[{),T], if^{s) ='*'(s), a.s. uj e il, 
s G [t,T], for some t G [0,r), then p(t; ^{■)) = "*"(•)), a.s. lo e Q. 

2) (Monotonicity) For any '!'(•), "*"(■) G L3^^[0,r], if ^{s) < "^(s), a.s. uj e VL, 
s G [t,r], for some t G [0,r), then p{.s;^{-)) > p(s;"^(-)), a.s. uj e n, s G [t,T]. 

Definition 2.3 A dynamic risk measure p : L'jr^[0,T] L|[0,T] is called a coherent 
risk measure if the following hold: 

1) There exists a deterministic integrable function r(-) such that for any ^{■) G 
L%[0,T], 

p(t;^(-) + c) =p(t;^'(-)) -ce^''^'"^"^'^", a.s. oj (£ ft, tG[0,r]. 

2) For ^i-) G L'^jr^[0,T] and A > 0, p{t;X^{-)) = Xp{t;^{-)) a.s. w G O, t G [0,r]. 

3) For any^{-),^i-)GL%JO,T], 

p{t;^{-)+^{-))<p{t;^{-)) + p{t;^{-)), a.s. w G 0, t e [0,T]. 
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Some necessary specifications on the generator g for BSVIE (2) are given by: 
(HI) Let g : x R"' x R""""^ x K"''"^ xQ he B{A^ x R"" x R"'""^ x R'^'"^)®J=-t- 

measurable sucli tliat s g{t, s, y, z, C) is J^-progressively measurable for all {t, y, z, () € 
[0,r] X R"" X i?'"^'^ X R""""^, and Vy, y € R"", z, z, C, C S R""""^, 

\9{t,s,y,z,C) - g{t,s,y,z,C)\ 
< Li{t, s)\y -y\+ L^it, s)\z -z\ + L^it, s)\C - C|, 
where {t,s) G A*^, Li{t,s) {i = 1,2,3) is deterministic non-negative functions. Further- 
more -E /q ( \gQ{t,s)\ds] dt < oo, where gQ{t, s) = g{t, s, 0,0,0). 

3 A new method for unique solvability of M-solution 

In this section, a new scheme is proposed and analyzed to simplify the unique solv- 
ability of M-solution in Yong [26]. The proof in [26] gives us a detailed procedure 
to comprehend how to construct M-solutions, however, from a mathematical point of 
view, it is rather tedious and sophisticated, and it should be of interest to develop a 
new brief approach for it. 

Inspired by the following equivalent norm for the elements of n^[0,T] in [22] . 



•H2[o,T] 



E r e%{t)\^dt + E r r e^'\z{t,s)\^dsdt 
Jo Jo Jo 



with /3 being a positive constant, we can propose a new one, see (6), and thus achieve 
the goal of giving a convenient and brief proof. In addition, compared with the proof 
in [22], it seems that the proof here is still simpler. Furthermore, we can also handle 
with the general case for p € (1,2] with this approach. 

Before doing this, some preparations are required. Consider the following simple 
BSVIE, 

Y{t) = -ip{t)+ [ h{t,s,Z{t,s))ds- f Z{t,s)dW{s). (4) 



(H2) h has the similar assumptions with g in (HI). Furthermore, L2{t,s) satisfies the 



condition, sup L2{t, s)^^''ds < oo, with some constant e > 0. 



te[o,T] 

The proof of the next proposition can be found in [26] . 



Proposition 3.1 Let (H2) hold, then for any € L^^[0,T], (4) admits a unique 
adapted M-solution {Y {■) , Z {■ , ■)) G nP[0,T]. Ifh also satisfies (H2), V^(-) G V^^[0,T], 
and {y{-), Z{-, •)) G T-L^[0,T] is the unique adapted M-solution of BSVIE (4) with (/i, ^) 
replaced by {h,ip), then Mt G [0,T], 

T \ f 1 



E { \Y{t) - Y{t)\P + ^ \Z{t, s) - Z{t, s)\''ds 
< CE |^(t) -■*(*) |P+ f / \h{t,s,Z{t,s)) -hit,s,Z{t,s))\ds 



(5) 
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Hereafter C is a generic positive constant which may be different from line to line. 
We move on to give the main result of this section. 
Theorem 3.1 Let (HI) hold, assume that 



sup 



t£[0,T] Jt 



T 



L'l{t,s)ds <oo, sup / L'^'^'^ {t, s)ds < oo, sup / L"?^ {t, s)ds < oo 



te[o,T] Jt 



where I + I = I, p e {1,2], ^ + ^ = 1, 1< p' < p. Then for any V'(-) € L^^[0,r], 
BSVIE (1) admits a unique adapted M-solution in W[{),T]. 

Proof. First let A^p[0,T] be the space of ah {y[-),z{-, ■)) € W[^,T] such that y{t) = 
Ey{t) + y^z{t,s)dW{s) with t G [0,r]. It is a matter of direct calculation to show 
that 7V1^[0, T] is a closed subspace of W[d, T] via the following two martingale moment 
inequalities in [5], 



E 



T 



z{t,s)dW{s) 



dt < CpE 



T 




* \ 2 

2 



z{t,s)\ ds dt, if p > 0, 



and 



E 







z{t,s)\^ds] dt<CpE 



z{t,s)dW{s) 



dt, if p > 1, 



where Cp is a constant depending on p. A new equivalent norm for the element in 
M'P[Q,T] of the form 



■)) II A1P[0, 



T] 



Jo 



E I e^'\y{t)\Pdt + E I e 



T 



\z{t,s)\^ds] dt 



(6) 



will be in force in the following part. Consider, 



Y{t) = iP{t) + / g{t, s, y{s),Z{t, s), z{s, t))ds 



Z{t,s)dW{s), 



(7) 



with t G [0,r], Vl-) G i^^[0,T] and {y{-),z{-,-)) G MP[Q,T]. Following Proposition 
3.1 we get that (7) admits a unique adapted M-solution {¥{■), Z{-, •)), and then define 
a map 6 : XP[0,r] ^ Mp[Q,T] by 

Q{y{-),z{;-)) = {Y{.),Z{.,-)), V(y(.),z(-,-))G-M^[0,r]. 

Let {y{-),z{-,-)) G MP[{),T] and e{y{-),z {■,■)) = (J {■) ,'Z {■,■)) , thus it follows from 
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inequality (5) that, 



E / e'''\Y{t)-Y{t)\Pdt + E / e 







p 

T _ \ f 

\Z{t,s) - Z{t,s)\^ds] dt 



< CE 



s, y{s), Z{t, s), z{s, t)) - g{t, s, y{s),Z{t, s), z{s, t))\ds \ dt 



< CE e^*|y^ Li{t,s)\y{s)-y{s)\ds^ dt 
+CE *{/ L3{t,s)\z{s,t) -z{s,t)\dsj dt 

P 

< CE [ e^* ( sup [ Ll{t,s)ds] [ \y{s) -yis)\Pdsdt 

Jo \t£\0,T]Jt I Jt 



+CE I e^* sup / Ll {t, s)ds 
\te[o,T] Jt J 



< CE \y{s) -y{s)\Pds e^^dt + C 



\z{s,t) -z{s,t)\P ds] dt 



p' 



T rT 



E I ds I e^'\z{s,t) -z{s,t)\Pdt 



< %E r e^'\y{s)-y{s)\Pds + C 



P Jo 

< ^E r eP'\y{s)-y{s)\Pds, 
P Jo 



p' 



^*dt / \z{t,s)-z{t,s)\Pds 



where l<p'<p, ^ + ^ = 1. Notice that here we use the following two relations, that 



p' ^ q' 

are, for any p G (1, 2], 1 < p' < p, r > 0, 

r-T 



\z{s,t) -z{s,t)\P ds 



< 



< 



' ds 



P-P 



P 



e'^'^\z{s,t)-z{s,t)\Pds 
^ e'"'^ [ e'^\zis,t)-zis,t)\Pds. 



(8) 



and E £ \z{t, s) - z{t, s)\Pds < CE\y{t) - y{t)\P which is a direct consequence of mar- 
tingale moment inequality and Holder inequality. Then we can choose a /3, so that the 
map is a contraction, and the result holds. □ 



4 A maximum principle for FBSVIE 

In this section, we give a stochastic maximum principle for forward-backward stochastic 
Volterra integral equations by assuming the control domain being convex and p = 2, 
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thereby generalizing for instance the case in [26]. As compared with the differential 
case, it is by no means clear that the method there can be extended to such setting. As 
we have claimed in the previous, there are some technical obstacles for us to overcome 
due to the absence of Ito formula here, in other words, we should adopt some other 
effective mathematical skills to circumvent the difficulties caused by it. Without loss 
of generality, we assume that m = d = 1. 



4.1 Setting the problem 

We denote by C/ a nonempty convex subset of R, and set 

U = {v{-) € L^O,T];v{t) £ U,a.s. t e [0,T],a.e.}. 

An element of U is called an admissible control. For any admissible control v{-) £ U, 
let us consider the following forward-backward stochastic Volterra integral equation, 
i.e., 

X{t) = ^{t)+ [ b{t,s,X{s),v{s))ds+ [ a{t,s,X{s),v{s))dW{s), 
< Jo Jo ^ 



Y{t) = m + / ait, s, y(s), Z{s, t),vis))ds 



T (9) 
Z{t,s)dW{s), 



associated with the cost functional by 



J{v{-))=E 



l{s, X{s),Y{s),v{s))ds + h{X{T)) + 7(y(0)) 







where ip{-) € Ljr[0,r] and ipi-) € Ljr^[0,T]. Basic assumptions imposed on 6, a, g, /, h, 
7 are stated as 



(H3) 



s, X, v) 
a(t, s, X, v) 
g{t,s,x,y,z,v) 
l{s,x,y,v) 



AxRxUxn^R, 
AxRxUxn^R, 
A^xRxRxRxUxn-^R, 
[0,T] X R X R X U X n ^ R, 



h{x) : X i? — > i?, 7(x) : Cl X R ^ R. b, a, g, I, h, and 7 are continuously differentiable 
with respect to the variables. The derivatives of b, a, g are bounded, the derivatives 
of / are bounded by C(l + \x\ + \y\ + \v\) and the derivatives of h, 7 with respect 
to X are bounded by C(l + \x\). Furthermore, we assume that gi{t, s,x,y, z,v) is 
;B(A'^x i?x i?x i?x C/)(8>J^T-measurable such that t gi{t, s, x, y, z, v) is J^-progressively 
measurable for all (s, x,y, z,v) £ [0,T] x R x R x R x U, {i = y,z). 

Given (H3) and v £ U, we observe that there exists a unique adapted M-solution 
{X{-),Y{-),Z {■,■)) G L3,[0,T] X L%[0,T] x ^^(0, T; L2,[0, T]) for above FBSVIE (9) 
by what we mean that X{-) satisfies the forward equation in (9) in the usual sense 
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and {Y{-), Z{-, •)) is the adapted M-solution of the backward form of (9). Both of the 
equations in (9) are called the state equations. The optimal control problem is to 
minimize the cost function J{v{-)) over admissible controls. An admissible control v{-) 
is called an optimal control if it attains the minimum. 

Remark 4.1 A special case of the above optimal control problem was considered in 
126^ where 



JK-)) = E 



l{s,X{s),v{s))ds + h{X{T)) 







and the coefficients are assumed to be independent of u. 



4.2 Variational equations and one convergence result 

Let u(-) be an optimal control and let {X{-),Y{-), Z{-, ■)) be the corresponding M- 
solution of (9). Let v{-) be such that u{-) + v{-) G U. Since U is convex, then for any 
< p < 1, Up = u{-) + pv{-) S U. Let's us consider, 

m = Mt) + [ bAt, s, X^{s),u{sms)ds + /* aS, s, X^{s),u{sms)dW{s), 







rj{t) = Mt) + gyit,s,X^{s),Y^is),Z''is,t),u{s)Us)ds 

+ [ g,it, s, X^{s),Y^{s), t),uis))as, t)ds - [ ({t, s)dWis) 
Jt Jt 



(10) 



where 



(^i(t)= [ b^{t,s,X''{s),u{s))v{s)ds+ [ a„{t,s,X''{s),u{s))v{s)dW{s), 
Jo^ Jo 

rT 



g^it, s, X^{s), Y^{s), Z^{s, t),u{sms)ds 
+ r g,{t, s, X^{s),Y^{s), Z^{s, t),u{s))v{s)ds. 



(11) 



The two equations in (10) are called variational equations. Obviously under assump- 
tion (H3) we can find a unique M-solution {^{■),v{-), C{-, ■)) G Ljr[0,T] x L'jr[0,T] x 
L^{0, T; Ljr[0, T]), which is the unique adapted M-solution of FBSVIE (10). We denote 
by {Xp{-),Yp{-), Zp{-, ■)) the M-solutions of (9) corresponding to Up. We now proceed 
to prove the relations 



e[ |Xp(t) -X"(t)pdt ^0; p^O, 

T 

e[ \Yp{t) -Y''{t)\'^dt ^0; p^O, 

T t 

E I |Zp(t,s) - Z"(t,s)pdsdt ^ 0; /) ^ 0. 
^ JO Jo 



(12) 
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In fact, it follows from the denotation of Xp, together with the forward equation in (9) 
that 

E / e-^'\Xp{t)-X'\t)\^dt 
Jo 

< CE [ e-^^dt [ \Xp{s)-X''{s)\^ds + CE f e'^^dt [ \up{s) - u{s)\'^ds 

Jo Jo Jo Jo 

< CE [ \Xp{s)- X''{s)\^ds f e-^^dt + CE [ \up{s) - u{s)\'^ds [ e"^^dt 

Jo Js Jo Js 

< -E [ e-^'\Xp{s)-X''{s)\^ds + -E [ e-^'\up{s) - u{s)\^ds, (13) 
1 Jo 1 Jo 



where 7 is a positive constant depending on the upper bound of all the derivatives. By 
choosing a 7 such that ^ = |, it leads to 

E r \Xp{t) - < e^'^E r e-^^\Xp{t) - ^0; p ^ 0. 

Jo Jo 

Following the conclusion of Theorem 3.7 in [26], we observe that 

e[ \Yp{t)-Y''{t)\^dt + E f [ \Zp{t,s) - Z''{t,s)\'^dsdt 
Jo Jo Jo 

< CE £ \g'it,s,Yp{s),Zp{s,t))-g''{t,s,Yp{s),Zp{s,t))\ds^ dt 

< CE \Xp{s) - X''{s)\'^ds + CE \up{s)-u{s)\'^ds^O, p ^ 0, (14) 

Jo Jo 

where g'{t, s, y, z) = g{t, s, Xp{s), y, z, Up{s)), g"{t, s, y, z) = g{t, s, X{s),y, z, u(s)), C is 
a constant depending on the upper bound of all the derivatives. Thus we can get (12). 
For t,s e [0,r], set 

Xp{t)=p-^{Xpit)-X-it))-m, 

Yp{t)=p-\Yp{t)-Y-{t))-rj{t), (15) 
_ Zp{s, t) = p~\Zp{s, t) - Z^{s, t)) - C{s, t). 

Using the similar method as (13), recalling the denotation of Xp, we can deduce that 
E r e-'''\Xp{t)\^dt<-E r e-'''\Xp{t)\'^dt + ep 

Jo OL Jq 



'P' 

where C is a constant depending on the upper bound of the derivatives, and Sp — ?> 0, 

QTlH 

a 2 ' 

rT 



yO — 7> 0. Then we can choose a such that — = i, and 



e[ \Xp{t)\^dt<e'^^E f e-"*|Xp(t)|2c/t < 2e"% ^ 0; /> ^ 0. 
Jo Jo 
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As to the term Yp, we arrive at 

E r e^'\Yp{t)\^dt + E r e^' f \Zp{t,s)\^dsdt 
Jo Jo Jt 

< CE [ e^* / \Xp{s)\^dsdt + CE f ef^f \Yp{s)\'^ dsdt 
Jo Jt Jo Jt 

+^E [ [ e^'\Zp{s,t)\''dsdt + Ce^'^e' 
P Jo Jt 

< ^E r eP^\Xp{s)\^ds + ^E r eP^\Yp{s)\'ds + CeP^ep, 
P Jo P Jo 



where C is an constant depending on the upper bound of ah the derivative, and £^—7-0, 
p — > 0. Then we can choose a /3 so that ^ < 1, and 

E [ e'^*\Yp{t)\'^dt + E [ e^* / \Zp{t, s)\'^dsdt 
Jo Jo Jt 

< CE [ e^'\Xp{s)fds + Ce^^e' 
Jo 

From above 

CE [ e^'\Xp{s)\^ds<Ce^'^E [ \X p{t)\'^ dt ^ 0; p ^ 0, 
Jo Jo 



thus ^ ^ 

E [ \Yp{t)\'^dt <E [ e^^\Yp{t)\'^dt ^0; p ^ 0. 
Jo Jo 

To sum up the argument above, we obtain: 

Lemma 4.1 Let (H3) hold, then 

f-T ^ fT _ 

limE / \Xp{s)\'^ds = 0, lim E / \Yp{s)\'^ds = 0. (16) 
p^o Jo Jo 

4.3 A simple form of stochastic maximum principle 

In what fohows, we make the following conventions with t,s € [0, T], f G U, 

l^is) = k{s,X'^is),Y^s),vis)),i = x,y,v, 

h^{s, t) = hj{s, t, X'" {t),v{t)), j = x,v, h = b, a, 

glis, t) = gk{s, t, X^{t), Y^{t), Z^{t, s),v{t)),k = x, y, z, v, 

where (X^, y, Z'") is the solution of (9) corespondent to v. In this subsection we assume 
that the cost function takes a simple form of J{u{-)) = E l{s, X{s),Y{s),u{s))ds. 
Since u is an optimal control, then p~^ [J{u + pv) — J{u)] > 0, and we have the following 
variational inequality. 
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Lemma 4.2 Let (H3) hold, then 



e[ l^{s)^{s)ds + E [ l''{s)r]{s)ds + E [ l'i{s)v{s)ds > 0, (17) 
Jo Jo Jo 

where {X, Y, Z) is the unique M-solution of FBSVIE (9) with u being an optimal control. 



Proof. From the Lemma 4.1, we know 



p-'E I [l{s,Xp{s),Yp{s),Upis))-l{s,X^{s),Y''{s),u{s))]ds 
E Us, X^is) + e{X,is) - X^is)),Y,is),Upis))^''^'^ ~ ^"^'^ 



ds 







P 



+E r lyis, X^{s),Y^{s) + eiYpis) - Y^{s)),up{s))^^^^^^^^^^ds 
Jo P 

+E I luis, X"(s), y"(s), u{s) + 0iup{s) - uis)))v{s)ds 
Jo 

E [ Us,X''{s),Y'\s),u{sms)ds + E [ ly{s, X^{s),Y'\s),u{s)Hs)ds 



+E / /„(s,X"(s),y"(s),'u(s))i;(s)ds. 
Jo 

Thus the conclusion follows. □ 

For deriving the maximum principle, we introduce the following two adjoint equa- 
tions: ^ ^ 

P{t) = l^{t)+ [ g'^{s,t)P{s)ds + [ g^As,t)P{s)dW{s), 

t T 

Q{t)=ll{t)+ f gl{s,t)P[s)ds+ f b%s,t)Qis)ds (18) 

+ / a^{s,t)R{s,t)ds - [ R{t,s)dW{s). 
Jt Jt 

Obviously the above FBSVIE admits a unique M-solution (P(-), Q(-), i?(-, •)) under 
assumption (H3). Note that gy and gz are non-anticipated processes under (H3). The 
later proposition is the so-called dual principle for linear stochastic Volterra integral 
equation, the proof of which can be found in 



Proposition 4.1 Let Ai : A x Q ^ R {i = 1,2) be B{A) (g) Tt -measurable such that 
s — )• A{t,s) is T -progressively measurable for all t € [0, T], furthermore, we assume 
that they are two bounded processes, ip{-) € Ljr[0,T] and € Ljr^[0,T]. Let ^(•) G 
L%[0,T] be the solution of FSVIE: 

i{t) = ^{t)+ f A^{t,s)i{s)ds+ f A2{t,s)i{s)dW{s), tG[0,r], 
Jo Jo 
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and {Y {■) , Z {■ , ■)) be the adapted M-solution to the following BSVIE, Vt G [0,r], 
Y{t)=i;{t)+ [ {Ai{s,t)Yis) + A2is,t)Z{s,t)}ds- [ Z{t,s)dWis). 



Then the following relation holds: 

rT 



E [ ^{t)i){t)dt = E [ ip{t)Y{t)dt. 
Jo Jo 



We now assert: 



Theorem 4.1 Let «(•) be an optimal control and y"(-), •)) be the corre- 

sponding M-solution of FBSVIE (9). Then we have, \/v € 

H{X''{t),Y''{t), Z'^it, ■),u{t),P{t), Q{t), R{-,t)) ■ {v - u{t)) > 0, a.e., a.s. 

where 

H{X-{t), Y^{t), Z-{t, ■),u{t), P{t), Q{t),R{; t)) 

= + QisK{s,t)ds + E^* [ R{s,t)a:^{s,t)ds+ [ g'i{s,t)P{s)ds 

Jt Jt Jo 

Here {P,Q,R) is the unique M-solution of FBSVIE (18). 

Proof. From the forward form in (10), the backward form in (18) and Proposition 4.1 
above, we know that 

E r mmdt+E r m f g:{s,t)p{s)dsdt 

Jo Jo Jo 

= E [ Q{t)dt [ b:^{t,s)v{s)ds + E [ Q{t)dt [ a^{t, s)v{s)dW{s) 
Jo Jo Jo Jo 



■■ E [ Q{t)dt [ 6^(t,s)u(s)ds + S / / R{t, s)a!iit, s)v{s)dsdt 
Jo Jo Jo Jo 

: E [ v{s)K{s)ds, (19) 
Jo 

K{s) = j\Q{tmt,s)+R{t,s)aUt,s)}dt. 

Similarly from the backward form in (10), the forward form in (18) and Proposition 
4.1, we know 

rT rT rT rT rT 

E / rj{t)q{t)dt = E P{t)dt / g^{t, s)Cis)ds + E / P{t)dt / g^{t, s)v{s)ds 
Jo Jo Jt Jo Jt 

rT rt rT rs 

= E ^{t)dt / gl{s,t)P{s)ds + E / v{s)ds / g^{t, s)P{t)dt. 
Jo Jo Jo Jo 

(20) 



where 
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It follows form (19) and (20) 

rT rT 



E f Cm{t)dt + E [ ii{t)l^y{t)dt 
Jo Jo 

pT r pT pT ps 

E / v{s)ds / Q{t)bl{t,s)dt+ / R{t,s)al{t,s)dt+ / g'!^{t,s)P{t)dt 

Jo Us Js Jo 



From the variational inequality (17) we have 

rT rT 



< e[ l^{t)^{t)dt + E [ l'^{t)ri{t)dt + E f l'^{t)v{t)dt 
Jo Jo Jo 

= E v{t)L{t)dt, 
Jo 

where 

L{t) = l'i{t)+ r Q{sK{s,t)ds+ r R{s,t)a'i{s,t)ds+ f g'i{s,t)P{s)ds. 
Jt Jt Jo 

The proof is complete. □ 

4.4 A general stochastic maximum principle 

In this subsection we will give a new maximum principle, while the cost function is a 
more general form 



J{v{-)) = E 



T 



l{s, Xis), Y{s),v{s))ds + h{X{T)) + 7(y(0)) 



It can be easily checked that E /"(s)^(s)(is + E Jq ly{s)r]{s)ds + E l'!^{s)v{s)ds + 
S/i^(X"(T))^(r)+£;73/(y"(0))?7(0) > O. in fact, the definition of implies E\Xp{T)\^ < 
6i{p) + CE \Xp{s)\'^ds, with 5i{p) 0, p —> 0, C is a constant depending on the 
upper bound of all the derivatives. Recalling the result in Lemma 4.1 we obtain 
E\Xp{T)\^ ^0,p^0. Similarly by the form of Yp, it follows that 

E\Yp{0)\'^ <CE [ \Xp{s)\^ds + CE f \Yp{s)\^ ds + C E f \Zp{s,Q)\^ds + 52{p), 
Jo Jo Jo 

with 62{p) — )■ 0, p — 7- 0. By the definition of M-solution in the previous, it does no 
matter what value of C('5,0) is as long as it is a deterministic function on s € [0,T]. In 
particular, if E \Dsr]{t)\'^ds < oo (here D is a malliavin operator, see jllj for more 
detailed accounts for malliavin calculus), then by Ocone-Clark formula (see [11]) and 
the definition of M-solution, we obtain 

r]{t) = ET]{t) + [ C{t,s)dW{s) = Er,{t)+ [ E^^ DsT]{t)dW{s), 
Jo Jo 
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thus we have C(^) -s) = E^^" Dsr]{t), then without loss of generahty we can determine 
C(s,0) by Er]{s). Similarly Zp{s,0) = EYp{s), and this leads to E\Yp{0)\'^ 0, with 

Summing up, we finally obtain p-^E{h{Xp{T)) - /i(X"(r))) ^ Eh^{X''{T))^{T), 
and p-^E{-f{Yp{0)) - 7(y"(0))) ^ E^y{Y''{0))r]{0), with E\C{T)\'^ < oo, E\t]{0)\'^ < 
oo, which is easy to validate. 

It follows from the martingale representation theorem that there exists a unique 
process Tr{s) G Lj,[0,T] so that /i^.(X"(r)) = ^/i^.(X"(r)) + Tr{s)dW{s); then 



Eh.iX^iTmT) 
EK{X^{T)) 

+Eh^{X%T)) 



Vi{T,s)v{s)ds+ / a'{;{T,s)v{s)dW{s) 



f bt{T,s)as)ds+ [ aUT,s)as)dW{s 
Jo 

= e[ b'i{T,s)h^{X''{T))v{s)ds + E [ 7r{s)a^,{T,s)v{s)ds 
Jo Jo 

+E [ bl{T,s)i{s)hAX''{T))ds + E f a^iT,s)as)7ris)ds. 
Jo Jo 

On the other hand, using the fact that Ejy{Y'^{0)) ({0, s)dW{s) = 0, one gets 

E7,(l^"(0))r?(0) = E r g:iO,shyiY^mas)ds + E g:iO,shyiY^mvis)ds 

Jo Jo 

+ e[ g^{0,shy{Y^{0)Ms)ds + E [ g]^{0, s)jy{Y^mEr]{s)ds. 
Jo Jo 

(21) 

In this case, FBSVIE (18) is replaced by 
f Pit) = l^it) +g-{0,t)^y{Y-iO)) +py{Y-iO))Eg-iO,t) 

+ [ g^{s,t)P{s)ds+ [ g^{s,t)P{s)dW{s), 
Jo Jo ^ 

Q{t) = lUt) + b^iT,t)h,iX^iT))+a^iT,t)7Tit)+gU0,thy{Y^m + [ g^is,t)Pis)ds 



+ j b^{s,t)Q{s)ds + j a^{s,t)R{s,t)ds - j R{t,s)dW{s). 

(22) 

So by a similar proof as Theorem 4.1 we get a more general stochastic maximum 
principle. 

Theorem 4.2 Letu{-) be an optimal control and {X{-),Y{-),Z {■,■)) be the correspond- 
ing M-solution of FBSVIE (9). Then we have, Mv S U, 

H{t, X"(t), y"(t), Z"(t, ■),u{t),P{t),Q{t), R{-,t)) ■ {v - u{t)) > 0, a.e.,a.s. (23) 
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where 

H{t, X{t),Y{t),Z{t, ■),u{t),P{t), Qit), R{;t)) 
= I'iit) + b':;{T,t)E^^h,{X^{T)) + a:iT,t)7T{t) + g:{0,t) 

+ / g'^{s,t)P{s)ds + E^' [ R{s,t)a'i{s,t)ds + E^' [ Q{s)b'^,{s,t)ds 
Jo Jt Jt 

where {P,Q,R) is the unique M-solution of FBSVIE (22). 
If we define 

-Hit, X^{t),Y^{t), ■),u{t),P{t), Qit),R{;t),v) = -H-v, 

then (23) can be rewritten as 

-Hit, X^{t),Y-{t), Z^{t, •), Q(i), R{;t), u{t)) 

= maxn{t,X''{t),Y''it),Z''{t,-),u{t),P{t),Q{t),R{;t),v) (24) 

We call Ti the Hamiltonian of the optimal control problem of FBSVIEs, and (24) the 
maximum principle condition. 

We would like to conclude this section by giving a application, that is, a linear- 
quadratic (LQ for short) problem of BSVIEs. The linear BSVIE is of the form 

Y{t)=i){t)+f [h{t,s)Y{s) + l2{t,s)v{s) + k{t,s)Z{s,t)]ds- [ Z{t,s)dWis), (25) 
Jt Jt 

while the cost functional J(^{-),v{-)) associated with the terminal condition ip{-) and 
control v{-) is given by 

rp 

J{i,{-),v{-)) = ]^E [Q'{t)YHt) + li{t)v\t)]dt+\EG'Y\Q). 

The linear-quadratic control problem is to minimize the cost function over admissible 
controls. Necessary assumptions will be in force in the following. 

(H4) Let /j : A*^ — )• R, {i = 1,2,3) be three continuous bounded processes such 
that s li{t,s) is F adapted for t € [0,T]. Q' and R are bounded and non-negative 
adapted processes, moreover, R'{t) > 5, where 5 is a positive constant, G' is a non- 
negative bounded random variable, Tp{-) € Lj-^[0, T]. In addition, assume U is also 
closed. 

Obviously (H4) is sufficient for the finiteness of the above linear-quadratic problem. 
Following the idea of Theorem 5.2 in Chapter 2 of [28], we are ready to present a 
existence theorem of the optimal control. 

Lemma 4.3 Let (H4) hold, then there exists au{-) € U such that J{ijj,u) = inf J{7p,v). 
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Proof. Let ^p{■) be fixed, and Uj{-) S -L|.[0,T] be a minimizing sequence of LQ problem, 
that is 

lim J(^(.),n,(-))= inf J(V(-),0). (26) 

Let {Y^,Z^) be the state processes corresponding to Uj{-). It follows from (26) that 
there exists a constant M such that J{ip{-),Uj{-)) < M for any j > 1. Additionally, 
J{ip{-),Uj{-)) > 6E Jg \uj{t)\'^dt, so we have E \uj{t)\'^dt < Consequently, there 
is a subsequence, which is still labeled by Uj{-), such that, 

Uj{-) u'{-), weakly in L|[0,T], 

By Mazur's theorem, we have a sequence of convex combinations 

Uj{-) = '^aijUi+j{-), Oij > = 1, 

i>l i>l 

such that 

Uj{-) —7- n'(-), strongly in L|[0, T]. 

Since the set U is convex and closed, it follows that u'{-) (zU. On the other hand, the 
Theorem 3.7 in 1261 leads to 



El \Yj{t)-Y'{t)\^dt<CE I / {uj{s) -u'{s))ds dt<CE \uj{s)-u'{s)\'^ds 
Jo Jo Ut J Jo 

i.e., Yj{ ) — )■ Y'{-), strongly in L|[0, T]. By the convexity of the generator for (25), 
Jm-),u'i-)) = hm J(^(.),%(-)) 



< lim S2aijJ{tp{-),Ui+j{-)) = inf J{ip{-),u{-)), 

which means that u'{-) is an optimal control. □ 

In this setting, the maximum principle condition can be written as 

rt 

.2/ 



-R!{t)u^{t) -l2{<d,t)u{t) -u{t) / l2{s,t)P{s)ds 

Jo 

> -R'{t)u{t)v -l2{0,t)v -V I l2{s,t)P{s)ds, (27) 

Jo 

with V e U, and this leads to R'(t)u{t) + ^2(0,^ + Jq l2{s,t)P{s)ds = 0, i.e., u{t) = 
R'~^{t)[l2{0,t) + f^l2{s,t)P{s)ds], where 

P{t) = Q'{t) + [h{0,t) + h{0,t)]G'Y{0)+ [ h{s,t)P{s)ds + [ k{s,t)P{s)dW{s). (28) 

Jo Jo 

Hence u{t) is the only control which satisfies the necessary conditions of optimality. By 
Lemma 4.3, it must be the unique optimal control. So we have 
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Theorem 4.3 Let (H4) hold, there is a unique optimal control n(-) for the linear- 
quadratic control problem. Moreover, u has the representation: u(t) = R'^^ (t)[l2{0,t) + 
Jq h{s,t)P{s)ds], where P{s) satisfies (28). 

By the form of the optimal control, we deduce that the optimal control indeed a 
linear state feedback of the entire past history of the state process P(-) instead of being 
a feedback of the current state, which is similar to the result for linear-quadratic control 
of BSDEs, see p.6-p.7 in [6]. Substituting the representation of n(-) into (25), together 
with equation (28), we get the following coupled FBSVIE 

P(t) = Q'(t) + [Zi (0,0 + ^3(0, t)]G'y(0)+ ! li{s,t)P{s)ds+ [ l3{s,t)P{s)dW{s), 
J, Jo ^ Jo 

y(t)=V(i)+ / l2is,t)R'^\s)l2iO,s)ds+ [ h{s,t)Y{s)ds 



T 



+ / l2{s,t) / l2{u,s)P{u)duds + / l3{s,t)Z{s,t)ds - / Z{t,s)dW{s), 
Jt Jo Jt Jt 

(29) 

Given (H4), by the unique existence of optimal control u(-), we observe that FBSVIE 
(29) admits a unique M-solution (y(-), Z(-, •), P(-)) by which means that P{-) solves 
the forward equation of (29) in the Ito sense and {Y{-), Z{-, •)) is the unique M-solution 
of the backward equation in (29). Consequently, 

Theorem 4.4 Let (H4) hold, then FBSVIE (29) admits a unique M-solution. 



5 Application in fincance 

In this section, we illustrate the maximum principle above by studying the risk min- 
imization problem in finance. Such kind of problem was studied by Mataramvura 
and 0ksendal [9] by formulating it as a zero-sum stochastic differential game. Re- 
cently 0ksendal and Sulem [12] also investigated this risk minimization problem via 
g-expectations. In this paper, we will consider the problem by means of the maximum 
principle in the previous. Some closed forms of the optimal solution are derived, which 
are consistent with the results in [12] or [23]. 

We consider a market with two investment possibilities, which are traded continu- 
ously until the fixed finite horizon T, is reached. One investment is described by 

dSo{t) = So{t)p{t)dt, So{0) = so. 

The other financial instrument is described by 

dSiit) = Siit)a{t) + Si{t)^{t)dW{t), Si{0) = si. 

Suppose the interest rate p{-) is nonnegative and bounded deterministic function, the 
stock-appreciation rate a(-) and the stock-volatility /3(-) are nonnegative and bounded 
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adapted processes. Moreover, /3 ^ {■) and {a{-) — p{-)) exist and bounded. The wealth 
process X{-) satisfies 

dX{t) = [p{t)X{t) + v{t){a{t) - p{t))]dt + v{t)P{t)dW{t), (30) 

with X{0) = X > 0, thereby the solution of the wealth equation can be given by 

X{t) = e^o Pi^y^x + / e-fs ''(")'^"[?;(s)(a(s) - + vis)P{s)dWis)]. 

Jo 

A portfolio v{-), representing the amount invested in the risk asset, is said to be ad- 
missible if v{-) € U. In the following the BSVIE that we are going to investigate is 

y(t) = -V'(t)+ / [ris)Y{s) + h{t,s)Z{s,t) + k2it,s)]ds- [ Z{t,s)dWis), (31) 
Jt Jt 

where ^ 

il^it) = h{X{T)) + ^ [h{t, s)X{s) + hit, s)v{s)]ds. 

Here r, li are bounded deterministic functions, ki is a process such that s — > ki{t,s) 
is F-progressively measurable for almost t E [0, T], ki is bounded, /c2 has the same 
assumption with go in (HI), h is utility function of the terminal wealth that satisfies 
(H3) and Eh?{X{T)) < oo (see the special cases below). If we define Q(t;ip{-)) = Y{t), 
then by Theorem 3.6 in Yong [25], g is a dynamic coherent risk measure. In order to 
show the result in a more explicit way, we would like to consider the special static case 
of g{0,ip{-)) = Y{0), and denote the cost functional by 

J{v{-)) = Y{0) = -Eij{0) + E [ [r{s)Yis) + ki{0,s)Z{s,0) + k2iO,s)]ds 

Jo 

= -Eh{X{T))-E [ [h{0,s)X{s) + l2{0,s)v{s)]ds 
Jo 

+E [ [r{s)Y{s) + ki{0,s)EY{s) + k2{0,s)]ds. 
Jo 

We want to find u & U such that J(n(-)) = inf J(v(-)). With the notation in the 

previous, we obtain 7 (y) = 0, and l{s, x,y,v) = [r{s) + Eki{0, s)]y — li{0, s)x — l2{0, s)v + 
k2{0, s). As to the coefficients in both (30) and (31), 

b:cit,s,x,v) = p{s),crx{t,s,x,v) = 0,K{t,s,x,v) = a{s) - p{s), 
a^{t,s,x,v) = P{s), g'^{t,s) = li{t,s), 

9'v{t,s) = l2{t,s),g'y{t,s) =r{s),g'^{t,s) = ki{t,s), 

where 

9i(t, s) = gi{t, s, X, y, z, v), i = x,v, y, z. 
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Then the Hamilton function is the form of 



n{t, y"(t), ■),u{t), p{t), Q{t),R{;t),v) 

-V [{a{t) - p{t))E^^h,{X^{T)) + m^t)] 



I l2{s,t)P{s)ds + p{t)E^' [ R{s,t)ds + {a{t) - p{t))E^' [ Q{s)ds 
'o Jt Jt . 

and the adjoint equation 

P{t) = [r{t) + Eki{0,t)] + r{t) [ P{s)ds+ [ ki{s,t)P{s)dW{s), 

Jo Jo 



Q{t) = -h{0,t) - p{t)h,{X^{T)) + / h{s,t)P{s)ds 

T T ° 

+ / p{t)Q{s)ds- [ R{t,s)dW{s), 



(32) 



h,{X^{T)) = Eh,{X^{T)) + / T,{s)dW{s). 

Jo 

Since is a hnear function in then the coefficient of v vanishes, i.e., 

= {a{t) - p{t))E^'h^{X''{T)) + /3(t)7r(t) + M{t), (33) 

where 

M{t) = [ l2{s,t)P{s)ds + P{t)E^' [ R{s,t)ds + {a{t) - p{t))E^' f Q{s)ds. 
Jo Jt Jt 

(34) 



The lemma below is essentially similar to Theorem A. 2 in the appendix of [12]. For 
readers' convenience, we present a proof here. 

Lemma 5.1 Consider the equation, a\{t)E-^*$^ + I3\{t)0{t) = C,{t), where ^ is a Tt- 
measurable random variable satisfying = E^ + 9{s)dW{s). Assume /3j~^(t) exists, 
C(t) is a adapted process, then ^ must be the form of 

^ = e-^^^'^^Ei + e-^i(^) r e^^^'^p^\s)as)dW{s), 

Jo 

where A\ is given by 

Mt) = f Pl\s)ai{s)dW{s) + \ f I3^\s)al{s)ds. (35) 
Jo ^ Jo 

Proof. We denote by P{t) = E^^i, therefore P{t) = P(0) + jl^e{s)dW{s). On the 
other hand, 

e{t) = f5^\t)m - a^{t)P{t)], 
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so ^ 

P(t) = P(0)+ / f3^\s)[Cis)-ai{s)P{s)]dW{s). 
Jo 

Since we can rewrite P by 

JO 

with Ai given by (35), thereby 

JO 

The conclusion foUows clearly. □ 

Remark 5.1 There are two things worthy to point out. On the one hand, if we assume 
that hx{X'''{T)) G Di,2, see JTl!, then by O cone- Clark formula, 7r{t) = E^Wth^iX"" {T)) , 
then (33) can he rewritten as 

{a{t) - p{t))E^^h.,{X'\T)) + mE^'Dth.AX^{T)) + M{t) = 0. 

It is a linear inhomogeneous Malliavin- differential type equation in the unknown random 
variable hx{X^{T)), which can also seen in llOij and HE/- On the other hand, if C = 0, 
there are infinite random variables satisfying the equation in Lemma 5.1. For example, 
if satisfies it, so does c^i, with c being a constant. 

To sum up, we give 

Theorem 5.1 Suppose u is an optimal portfolio of the above risk minimizing problem, 
then u{-) must satisfies 

i-T 
Jo 

with X"(r) be the terminal wealth corresponding to u, M{t) given by (34), and 

A{t)= f f]-\s){a{s)-p{s))dW{s) + }- f !3-\s){a{s)-p{s)fds. 
Jo ^ Jo 

In order to express the explicit form of u, some more assumptions are required. Let 
li = 0, (i = 1, 2), then the Hamilton function becomes 

U = -v[{a{t)- p{t))E^^h,{X'\T)) + m<t)\ 

-V (^p{t)E^' R{s,t)ds + {a{t) - p{t))E^' Q{s)ds^ , 
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where 



Q{t) = -p{t)h^iX''{T))+ [ p{t)Q{s)ds- [ R{t,s)dW{s) 



h,{X^{T)) = Eh^{X^{T)) + [ 7r{s)dW{s), 

Jo 

therefore, the optimal portfoho u satisfies 

= {ait)-pit))E^'h,iX^{T))+mAt) 

+/3{t)E^' R{s, t)ds + (a(t) - p{t))E^' Q{s)ds. (36) 

As a consequence, by solving the above simple BSVIE, we deduce that V(t, s) E A, i.e., 
< s < t < T, 

Q{t) = -p{t)e^^ P^^)'^^ ■ E^*h^{X''{T)),R{t,s) = -p(t)e/t^ ''(")'^"7r(s). (37) 

Substituting (37) into (36), we arrive at 

{a{t) - p{t))E^^K{X{T)) + m<t) = 0. (38) 

Recalling Lemma 5.1, we can express hx{X^{T)) as /ia;(X"(r)) = e'^^'^'^'i EK{X''{T)), 
which is a necessary condition for u being optimal. 

Remark 5.2 On the one hand, due to (38) being a homogeneous Malliavin- differential 
type equation, ifh{-) is replaced with ch{-) in the cost functional, with c being a constant, 
we can still obtain the same result. On the other hand, thanks to I2 = 0, the Hamilton 
function % is independent of P, which is solution of the forward adjoint equation in 
(32). Then no matter what values ofki, it does not change the value of optimal portfolio 
u. 

Now we will prove that the necessary condition above is also sufficient. For any 
Vi^U with i = 1, 2, we have form the concavity of h that Eh{X''^ (T)) - Eh{X''^ (T)) > 
£^[/ia;(X^i(r))(X''i(r) - X''2(r))], where X''' (r) is the terminal wealth corresponding 
to Vi, thus one sufficient condition for the strategy u being optimal is that E[hx{X'^ {TYiX"" {T)] 
being a constant over v (zU, see Proposition 2.1 in Z. Wang [23]. By the above neces- 
sary condition we have Ehx{X''{T))X''{T) = £;/i^(X"(r)) • Ee-'^^'^'^X"{T). Using Ito 

formula to e-^Wx^(^)ei't'^''(")'^^ one gets 



A(s)g/Jp{s)a 



v{s)m-X{s)^^ 



dW{s), 



hence Ee''^^'^') X" {T) = xe^o^ £;/i^(X"(T))X^(T) is a constant independent of v. 
Thus we obtain 
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Theorem 5.2 Suppose u is an optimal portfolio of the above risk minimizing problem 
if and only if u{-) satisfies h^{X'^{T)) = e-^(^)£;/i^(X"(r)). 



Remark 5.3 There is one thing worthy to point out. The above argument also implies 
that the sufficient condition in Proposition 2.1 in f2^ is also necessary. In fact, when 
h is a concave utility function, from our stochastic maximum principle we know that 
one necessary condition for u being optimal is equation (38), and this also implies that 
Ee~^^'^^ X"" {T) is a constant independent of v. 

Some special cases for the function h{x) are given below to show some exact results 
of u. 

Case 1 h{x) = x, we will deduce that the optimal portfolio is u{t) = 0. In fact, 
from equation (38), we have a{t) — p{t) = 0, then a{t) = p{t), which means the stock- 
appreciation rate is equal to the interest rate. In this case, the optimal portfolio is 
u{t) = 0. As we know, there are usually many kinds of method, i.e., risk measure, to 
measure the wealth at some time, for example the terminal wealth at time T, and there 
is one optimal portfolio for each kind of risk measure. On the other hand, from Remark 
5.2, we can choose any bounded ki{t, s), in other words different risk measures, thus get 
different value Y{0), i.e., different minimal risk, while the optimal portfolio is the same 
one. For example, if ki{t, s) = 0, k2{t, s) is independent of t, then the minimal risk can 
be expressed as, Y{0) = e/o^/'(^)+^(^)°'^x + E /g^ e/o K«)'i«fc2(s)ds. Note that p = 0, it is 
consistent with the one in 



Remark 5.4 Recently, the author ]20^ consider the case when r is allowed to be ran- 
dom, while r is assumed to be deterministic in 125^ . In this case, it is easy to show that 
the above results also hold, and the minimal risk is given by y(0) = Ee-I'o P(^)+K^)'i^ 

X + 

E Jq el'o^^^)'i^k2{s)ds. 

Case 2 h{x) = x — with 7 7^ being a constant, then hx{x) = 1 — jx, and in the 
following, we denote F{t) = E^^hxiX'^iT)) = l--fE^^X{T). In this setting, we assume 
thatE Jo \v{s)\''ds < 00. Following the idea in [23], we will show the explicit form of the 
optimal portfolio u. By equation (38), we know that 7r(t) = —f5"^{t)F{t)[a{t) — p{t)). 
Using Ito formula to A{t)F{t) on [0,T], where A{t) = e^o '^(«)'^« and a is a deterministic 
integral function, 

A{T)F{T) = F(0) + f A{s)T:{s)dW{s) + f F{s)a{s)A{s)ds. (39) 
Jo Jo 

Since F{T) = 1 - 7X"(r), together with (39), we have 
X^{T) = - (^A(T))"i ' 



F(0)+ [ A{s)TT{s)dW{s) + [ F{s)a{s)A{s)ds 
Jo Jo 



-< -^)-L iA(r) "'^L — ^Mfm — 
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On the other hand, by (30), 



T 



X"(T) = e/o Pi')'i'x + / e-fs ^'(")'^"(q(s) - p{s))u{s)ds + ''(")'^"/3(s)u(s)dTy(s) ,(40) 



'- 

then by comparing the correspondent part in (39) and (40), 



p{s)ds 



F(s)a(s)A(s) _ 
7A(T) - 

A{s)F(s){c{s)-p{s)) 
lA{T)l3(s) 



p{u)du 

{a{s) - p{s))u{s), 



T 



p{u)du 
s f3{s)u{s), 



thereby we deduce that a{s) = , and the optimal portfoho is expressed as 

u{s) = 



- f p{u)du - p{s))F{s)A{s) 

P^{s)^A{T) 
-Sl O'''%^!:''" +pi-))dn {a{s)-p{s)) 

/32(s) 
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(41) 



To get more feehng about the general form of equation (33) in the previous, i.e., 
the inhomogeneous Malliavin differential equation, we will consider some special cases. 
Let li{0,t) = li{t), i = 1,2. Moreover, for the sake of obtaining the exact expression of 
P, Q, R, we assume ki{t,s) = r{s). 

Case 1 If p{t) = 0, then for any t € [0, T], s < t, 



Pit) 

Q{t) 

R{t,s) 



2r{t)e^o ^(*)'^*-|/o ''^(■5)'^«+/o r{s)dW{s)^ 
-hit) + hit) [ P{s)ds 



hit) (^-1 + 2 r(s)e/o '■(«)rf«-5/o '■'W'^^+Zo ^•H'^^Hfis 



then 



M{t) = 2l2it) / r(s)e^oK«)'i«-|/o r-2(«)d«+/Q=r-(u)dM/(n)^^ 

Jo 

+2Pit)E^' his) (^e^o - e/oK")'^"^ r(t)e-5 /o '■'{")'^"+/o '^W'^'^Wfis 
+ait)E^' r his) (-1 + 2 r riu)e^o''^'''^'^''-yo'^H^)'i^+foHv)dW{v)^\ 
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In fact, in this case, ly{s, x,y,v) = 2r(s), and P(-) and Q{-) satisfy 

P{t) = 2r{t) +r{t) [ P{s)ds + r{t) [ P{s)dW{s), 

< t T 

Q(t) = -li{t) + h{t) [ P{s)ds- [ R{t,s)dW{s). 
Jo Jt 

By the martingale representation theorem, there exists a unique adapted process such 
that 



P{t) = EP{t)+ [ T{t,s)dW{i 
Jo 



thus we have 



R{t,s) = hit) / T{u,s)du 



Assume r 7^ 0, and p'{t) = then we have 

P'{t) = 2+ [ r{s)P'{s)ds+ [ r{s)P'{s)dW{s), 
Jo Jo 

Using Ito formula to P"{t) = e" ^' and we obtain P"{t) = 2+ r{s)P"{s)dW{s), 

thus we solve that P"{t) = 2e-|/o ^'(^)'^^+/o ^(^)'^^(^), so we obtain P{t) above. On the 
other hand, EP"{t) = 2, then P"{t) = EP"{t) + Jlr{s)P"{s)dW{s), thus 

^- Slr{u)duPif) ^ ^-!lr{u)duEm ^ f\^s)P" {s)dW {s) , 

r{t) r{t) Jq 

thus 

T{t,s) = 2r(t)r(s)e-^o'"(")'^"-|/o ''^(")'^"+/o ''W'^^Wj 

thus we get 

R{t,s) = 2^i(t)r(s)e"^-^0 ^'W'^«+/o Lfor{v)dv _ ^Jg r{v)dv 

As to the case r = 0, then V(t, s) G A, P{t) = 0, Q{t) = -li{t), R{t,s) = 0, and they 
are all consistent with the above results. 

Case 2 If r(s) = 0, then we have P{t) = 0, and 

rT 



Q{t) = -h{t)-p{t)h,{X''{T))+ p{t)Q{s)ds- R{t,s)dW{s). 

Jt Jt 

and we can solve M-solution as follows 

Q{t) = -/5(i)e/*^''(")'^" •S^*/i^(X"(r)) 



25 



and 

R{t, s) = -p{t)e^t pW'^^7r(s), {t, s) G A, 

thus 

M{t) = (/3(t)7r(t) + (a(t)-p(i))i?-^*/i.(X"(T))) (l-e/*^''(^)'^^) 

-{a{t) - p{t))E^' (^h{s) + p{s) J\j's P^'''^'^''h{v)dv^ ds. 

Case 3 r(s) 7^ 0, p{s) ^ 0. In this case, we can obtain the following result by 
combining the results in above two cases together, 

^ 2^^-l^^Jor{s)ds-^ f^r'^{s)ds+f^r(s)dW(s) ^ 

Q{t) = Q'{t) - p{t)e^t Pi^)<i^ . E^'h^{X{T)) + p{t)E^' j\it p^'^^'^'^Q' {s)ds, 
R{t, s) = R'{t, s) - p{t)e^t ^(")'^«7r(s) + p{t) e^t p^-")'^^ Rl {u, s)du, 
where 

Q'{t) = h{t) (^-1 + 2 r{s)e^O r^du-^ r\u)du+J^ riu)dWiu)^^^ 

R'{t,s) = 2h{t) (^e^'^(")'^"-e^oK«)rf«^^(s)e-|/o'^'W<i«+/o'^W'iWW, 

thus 

M{t) = {a{t) - p{t))E^'h^{X''{T)) [l - eStPi^)d^^ + /3(t)7r(t) (l - e^t p(^^)d^^ 
+{a{t) - p{t))E^' {q\s) + p{s) e/>(")'^"Q'(?;)(i7;j ds 

+(5{t)E^' (^R'{s, t) + p{s) e^s Pi^)'^^R'{v, t)dv^ ds. 
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